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Viscous propagation of mass flow variabiflty in accretion discs 
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ABSTRACT 

We study mass flow rate through a disc resulting from a varying mass supply rate. Variable 
mass supply rate occurs, e.g., during disc state transitions, and in interacting eccentric bina- 
ries. It is, however, damped by the viscosity of the disc. Here, we calculate this damping in 
detail. We derive an analytical description of the propagation of the flow rate using the solu- 
tion of Lynden-Bell & Pringle, in which the disc is assumed to extend to infinity. In particular, 
we derive the accretion-rate Green's function, and its Fourier transform, which gives the frac- 
tional damping at a given variability frequency. We then compare this model to that of a finite 
disc with the mass supply at its outer edge. We find significant differences with respect to the 
infinite disc solution, which we find to overestimate the viscous damping. In particular, the 
asymptotic form of the Green's function is power-law for the infinite disc and exponential for 
the finite one. We then find a simple fitting form for the latter, and also calculate its Fourier 
transform. In general, the damping becomes very strong when the viscous time at the outer 
edge of the disc becomes longer than the modulation time scale. We apply our results to a 
number of astrophysical systems. We find the effect is much stronger in low-mass X-ray bina- 
ries, where the disc size is comparable to that of the Roche lobe, than in high-mass binaries, 
where the wind-fed disc can have a much smaller size. 

Key words: accretion, accretion discs - binaries: general - stars: individual: (4U 1820-303, 
Cir X-1, Cyg X-1 GX 301-2) - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 

If a binary in which mass is transferred from one component to 
the other is ecce ntric, the mass transfer rate will be periodic at 
the bi n ary period i Avni. F abian. & Pringle 1976; Avni & Goldman 
198Cl; iHavnes. Lerche & Wright 1980.; Hut & Paczvnski. 1984 ; 



Brown (feBoyld Il984l; lBovle& Walked Il986l; IStevenj 1 19881 ; 



Lavton et aljil998l ; iRegos. Bailey & Mardling||2005l) . The larges t 
class of such bina r ies is Be/X- r ay bin aries, see, e.g., ICo3 l llOOOl) . 
IZiolkowskil ( l2002l) . lNegueruelal | |2004|) for reviews. All known sys- 
tems of this type consist of a (strongly magnetized) neutron star 
orbiting, usually in an eccentric orbit, a high-mass (unevolved) Be 
star. The accretion occurs via stellar wind and interaction with a 
circumstellar decretion disc a round the massive star, and is maxi- 
mized around the periastron dWaters et al.lll989 ). Indeed, most of 
those systems show X-ray p eriodicity at their orbital perio ds, so- 
called type-I X-ray outbursts jStella, White. & Rosneilll986l) . seen, 
e.g., in A0538-66, EXO 2030-1-375, A0535+26, AX J0049.4-7323, 
4U 0115+63 (see, e.g., references in Ziolkowski 2002). The accre- 
tion onto the neutron star in thos e systems is thought to occur via 
a (non- stationary) accretion disc jCorbelll 19861 ; IWaters et alJI 19891 ; 
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iHavasaki & Okazakill2004l) . Arguments for the existence of an ac- 
cretion disc include a lag between the onset of optical activity and 
an X-ray outburst, abrupt rises and linear decays of outburst (im- 
plying the storage of materi al in a reservo ir), and spin-up of the 
neutron star during outburst ( ICorbel|[l986l) . The last phenomenon 
is seen during X-ray pulsations at the neutron star spin frequency. 

Another class of binaries relevant to this work consists of X- 
ray pulsars with supergiant donors. Here, accretion is via stellar 
winds and, in some cases, via a Roche lobe overflow. Some of those 
systems have eccentric orbits, which then cause p eriodically vary- 
ing accretion, e.g. , GX 3 01-2 H Kelley. Rappap ort. & Petrdll98d; 
iMukheriee & Paull |2004| ; iLeahv & Kostka ,20081) . It has the or- 
bital period of Pprb - 41.5 d and t he eccentricity of e 0.5 
jKaper. van der Meer & Naiarroll200^ . The presence of a n accre- 
tion disc in this system is discussed by IPravdo & GhoshI ( 1200 ih . 
and details of the mass transfer from the compani on, including the 
presence of an accretion stream, are studied by iLeahv & Kostkal 
(2008). 

Another eccentric X-ray binary is the peculiar system 
Cir X-1, which contains either a low or high-mass donor 
djonker. Nelemans. & Bass3 20071). of Pprb - 16.5 d and e ^ 
0.5 ('Parki nson et alj |2003| ; ]wen et all |2006| ; Ijonker et al.l |2007|) . 
The sy stem shows strong periodicity at the orbital period in 
X-rays dParkinson et"aLll2003l ; Iciarkson. Charles. & Onvet3l2004l; 



© 2009 RAS 



2 A. A. Zdziarski, R. Kawabata and S. Mineshige 



IWenetal.|[2003) and in radio dWhelanetalllig?-?!) . Accretion in 
that system may occur vi a a Roche-lobe overflow near the perias- 
tron dClarkson et alj2004t) . 

Then there is the case of the ultracompact low-mass X-ray 
binary 4U 1820-303, where P^rb - 685 s is detected in X-rays 
jStella. Priedhorsky & Whi"t3 ll987: Zdzia rski et alJ l2007b^. In ad- 
dition to the bina ry period, it shows a strong periodicity at 170 d 
dWen et alj|2003) , and there is compelling evidence that it is due 
to the corresponding long-ter m modulation of the accretion rate 
dChou & Grindlavl l2001 : Zdzi arski. Wen & Gierhhskill2007iJ) . The 
main model of the long-ter m modulation is an int eraction with a 
third star ( iChou & Grindiavll200ll : IZdziarski et al.ir2007a) . which 
causes a peri odic modulation of the eccentricity of the inner binary 
( lKozailll962h . Though e <K 1 in this model, the mass supply rate to 
the outer disc is strongly modulated at both the orbital and superor- 
bital periods. An important issue is whether thi s orbital modulation 
of the mass supply rate can b e seen in X-rays dOilfanov & ArefievI 
l2005l : IZdziarski et al.ll2007bh . 

Enhancement of the mass-transfer rate around the perias- 
tron in low- mass X-ray binaries with e <k 1 has also been 
proposed by iMaccaroM l l2005h to be responsible for flares 
discovered from some of X-ray sources in NGC 7469 by 
ISivakoff. Sarazin. & JordanI j2005h . Again, an issue to examine 
is whether periodic enhancements in the mass transfer rate in 
those sources can actually lead to X-ray flares. Viscous damp- 
ing of a change in the mass supply rate is also an important is- 
sue for understanding stat e transitions in X-ray binaries, see, e.g., 
iHameurv & Lasot j ll2002f) for the case of a cataclysmic variable. 

Thus, the question of the response of an inner part of the disc 
to a changing mass supply rate in its outer parts is an important 
astrophysical issue. Here, we study this problem by calculating the 
effect of the disc viscosity on the propagation of a disturbance from 
the disc outer edge to the inner one. Our main method (Section[2} 
is to calculate Green's functions for the problem, i.e., the disc re- 
sponse to a i5-function pulse of the mass supply. Convolution of the 
Green's functions with any form of mass supply rate, see equation 
([8}, gives the mass accretion rate at the inner edge of the disc. The 
squared Fourier transform of the Green's function multiplied by the 
power spectrum spectrum of the mass supply rate gives the power 
spectrum of the accretion rate, see equation Q. In particular, if the 
supply rate is sinusoidal, the absolute value of the Fourier transform 
gives the factor by which the original amplitude is reduced. After 
presenting our theoretical results (Sections |2l [3}, we apply them to 
X-ray binaries (Section|4l(. 



2 THE GREEN'S FUNCTIONS FOR THE ACCRETION 
RATE AND THEIR FOURIER TRANSFORMS 

We calculate propagation of variability of the mass flow rate 
through an accretion discs by solving time-dependent disc struc- 
ture. We assume that the disc is geometrically thin, axisymmet- 
ric and Keplerian. The equations o f continuity and angular mo- 
mentum transport by viscosity are jLvnden-Bell & Pringldri974 
iKato, Fukue & Mineshig3ll998h . 



Ik 
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1 dM 

Inr dr 
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(1) 
(2) 



where r, E and M are the disc radial coordinate, the surface density, 
and the mass flow rate, respectively. Instead of solving the energy 



equation, we assume that the kinematic viscosity, v, is a power-law 
function of only r. 



V = v(rou,) — , 
\ '•out / 



(3) 



where Tom is the disc outer radius. In the g- viscosity prescrip - 
tion, y = (2/3)c!'cJ/f2K (using the conventions of lKato et alJI 1998b . 
where a is the viscosity parameter, Cs is the sound speed, and 
is the Keplerian angular velocity. Then, possible values of n in- 
clude 1/2 for a disc without cooling, 3/2 for an isothermal disc, 
and 3/5, 3/4 for a thin disc supported by gas pressure and domi- 
nated by eith er electron scattering or bo und-free Kramers' opacity, 
respectively jShakura & Sunyaevlll973h . In general, the viscosity 
as a function of radius may not be a power law, e.g., when the dom- 
inant source of opacity changes with the radius. We also note that a 
propagating disturbance in S will change the local temperature, and 
thus both th e sound speed and visco s ity, v. How ever, following th e 
treatment of iLynden-BeU & Pringl3 ( 1 19741) and lKato et al] ( Il998h . 
we neglect here this nonlinearity. 



2.1 The analytical model of an infinite disc 

The propagation of variability can be calculated by a simple ana- 
lytical model assuming the disc extending from infinity to r = 0, 
where the density vanishes. At an initial moment, f = 0, the col- 
umn density is given by a delta function, X = So<5(^ - 1), where 
if = (r/ro)''^. Here, ro is the radius at which mass is supplied to the 
disc. We will compare in Section l2!2l below the analytical results 
assuming the disc extending to infinity and the mass injection hap- 
pening at a finite radius to the numerical ones, where we assume 
the mass is supplied at the outer edge, Tout, of a finite disc. With 
the above assumption for the mass injection, the solution of equa- 
tions l[l}{3} constitutes th e Green's functions for the problem. It can 
be ob tained analytically jLynden- BeU&Pringl3ll974 l Kato et al.l 
1 19981) . 



+ 1) 
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where Im(z) is the modified Bessel function of the first kind. 



(4) 
(5) 



1 



A-2n 



2r^ 

(6) 



the viscous time scale is defined as fyisc = rl\vi\ (for a steady-state 
disc), and v,- is the radial velocity. Hereafter, we will consider the 
cases with n < 2, corresponding to n > Q, tho ugh the solutions 
for n > 2 can be obtained as well (Kato et al.lll99S) . At n < 2, 
?visc increases with the increasing r, see equation {6]l. The Green's 
function for the c olumn density propagation ha s been used, e.g., b y 
iMineshige & Wo od ( 1989), Lyubarskii ( 1997l), lKotov et alj ( l200lh . 
whereas that of equation (|5} appears not to have been used yet in 
the literature. 

Equation ([5j can be expressed in an explicit form at any radius, 
see Appendix |A] Here we are interested in the accretion rate at 
the disc inner edge. Hence, we consider the rate at ^ = 0, which 
Green's function we denote as M(,G^{t), where is defined to 
be dimensionless. We find. 



Gm(t) = -prr-r "exp 



(7) 
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Figure 1. (a) The Green's function for the profile of the accretion rate at 
the inner edge of a disc for an infinite disc with the mass supply at a finite 
radius, in the cases of n = 3/4 (red solid curve) and n = 1 (blue dashes). The 
black dots (resolved only at small values of r) show the numerical Green's 
function for the finite disc at n = 1 , see Section [T2l (b) The absolute values 
of the (infinite-disc) Fourier transforms for n = 3/5 (black dots), n = 3/4 
(red solid curve) and « = 1 (blue dashes). 



where T is the Euler Gamma function, GM{T)dT = 1, and 
G^j(t < 0) = 0. This rate is plotted in Fig. [fja) for n = 3/4 
= 2/5) and « = 1 = 1/2). 
If the mass supply rate at rg is M„^t('^), the rate at the inner 
disc edge is a convolution, 



Mi„(r) : 



Moui(t - 

Jo 



Ti)GM(ri)dTi 



(8) 



The Fourier transform of /Wi„(?) is then the product of the individual 
transforms, i.e., 



(9) 



The Fourier transform, G]i](T)e "^'"dr, of the Green's function 
of equation ([7]l is. 



(11) 



where K„,(z) is the modified Bessel function of the second kind, 
ll) = 27r/?visc(?'o): and / is the frequency. For n = I (ji = 1/2), the 
transform assumes a particularly simple form, 

^MM = e-<'^"^, \gMico)\ = e-^. (12) 

Fig.IIJb) shows \ff^(a>)\. If the original signal is sinusoidal, its am- 
plitude is reduced by \ffM(<^)\- We see that the viscous damping in- 



creases with the increasing n. This is caused by the viscous time 
within the disc, t,,,,(r) = tyi^(r„M)(r/r„^tf'", being longer for a 
larger n at a given ?visc('"out)- 

In order to account for the effect of viscous diffusion on the 
power spectrum of an arbitrary mass supply rate, we need to mul- 
tiply its power spectrum by In general, damping (time- 
dependent or stationary) of a variable component of the accretion 
rate due to the viscous transport through the disc can be calculated 
using the convolution of the light curves, equation {8j. Damping 
can also be calculated by integration of resulting power spectra, 
equation Using the former method to calculate the time evolu- 
tion for a step-function mass-supply either increase from to Mq 
or decrease from Mq to at t = yields at r > 0, 



MJr) = Mof(T), M„,(T) = Mo[l - /(t)]. 



respectively, where 



fir)-- 



exp(-2/T), 



yu > 0; 

yU = 1 (« : 



3/2), 



(13) 



(14) 



and T(x,y) is the incomplete Gamma function. By subtraction, we 
can also calculate M-^^ for a rectangular pulse of Mout- Note that our 
results do not depend on the normalization of M as the equations 
are linear. In particular, a constant M can be added to any of the 
solutions. 



"~ 2.2 Finite disc 



In our numerical calculations, we assume the disc is cut off" at Tom, 
at which the mass accretion rate MoutCO is given as a boundary con- 
dition, and its density is zero at the inner edge, ri„. We also assume 
conservation of angular momentum at ro^, and do not allow a disc 
formation beyond it. Since the kinematic viscosity is assumed to 
be a power-law function of r, the results depend only on the ra- 
tio Tout/n,,. We calculate the time evolution of the disc structure by 
solving equations l[T}[3). We divide the disc radius into N = 10^- 
10' grid points. 

In order to calculate Green's function for this problem, we 
consider a pulse with the duration of 10"^fvisc(''out) at the outer edge 
of a disc with rout/^in = 10^, for n = 3/4, 3/5 and 1. We now define 
T = r/fvisc(''oui)- We have found out that the outer boundary con- 
dition changed with respect to that in the analytical model results 
in a rather significant change of the shape of the Green's function. 
Namely, whereas the analytical Green's function has a power-law 
asymptotic (at t » 1) form, Gj^(t) oc t"'"'', equation l|7]l, the nu- 
merical Green's function behaves like Gm(t) oc exp(-/3r), where 
we have found the dependence. 



1 1 
4 2^' 



(15) 



On the other hand, the shape (but not the normalization) of the 
finite-disc Green's function at low t is found to be virtually identi- 
cal to that of equation Combining these two dependencies, we 
find excellent fits to our numerical results by. 



, J^--exp(-^), 



T < To; 
r > To. 



(16) 



Here, the only free fitting parameter is tq. The value of rj follows 
from the continuity at tq, and A' follows from the integral over t be- 
ing unity. Both conditions are imposed during fitting (of InG^^). 
For completeness, we give here all the obtained values, tq = 0.781, 
0.400, 0.301, N = 2.190, 2.139, 2.296, t] = 0.901, 0,844, 0.669, 
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Figure 2. The Green's functions for the profile of the accretion rate at the 
inner edge of a finite disc with the mass supply at the outer edge. The black 
dots (resolved only at small values of r) give the numerical solution, and the 
green dots and the red solid curve give the fit of equation )16) for n = 3/5 
and 3/4, respectively. The results are compared to those for the infitite disc 
at n = 3/4, shown by the blue dashes. 
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Figure 3. Comparison of the absolute values of the Fourier transform and its 
approximation for n = 3/4. The black solid curve gives both the transform 
of the numerical Green's function for the finite disc and the transform of 
equation )16t , which are indistinguishable from each other on this plot. The 
green dots give the approximation of equation U7t . The blue dashes give 
the analytical result for the infinite disc, equation jlOt . 



for n = 3/5, 3/4, 1, respectively. We see that the initial (at low 
t) Green's function is larger by /V ~ 2 than that for the infinite 
disc. This can be explained by the propagation of the same M being 
now only inward. The Green's functions and their fits for n = 3/5 
and 3/4 are plotted in Fig. |2] We see that equation l ll6t provides 
excellent fits to the shown numerical results, within < 4 per cent 
at T>0.05. On the other hand, the shape of equation ^ gives a 
slightly worse description of the Green's function at t < To for 
n = 1, which value of n, however, does not correspond to any phys- 
ical disc. Figs.[TJa) and|2]compare the finite-disc Green's function 
to those for the infinite disc, equation l[7]l. 

The Fourier transform of the above Green's functions, see Fig. 
[3] can be obtained either directly from the numerical results or by 
integration of equation l |16t . The Green's function and its transform 
obtained here for a finite disc can be used analogously to those for 
the infinite disc. Section [211 e.g., to solve a time-dependent disc 
behaviour, damping of a periodic variability, or the efi'ect on power 
spectra. 

The narrowness of the Green's functions for the finite disc re- 
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Figure 4. The Green's function for the profile of the accretion rate at the 
inner edge of a finite disc with the mass supply at ro <K r„at for n = 3/4. 
The results are compared to those for an infinite disc, equation (7), shown 
by the thick solid curve. 



suits in the transform for the finite disc being larger by ~ 2 at high 
/ than that for the infinite disc, though both equal 1 at / — » 0, see 
Fig. [3] We have found that the ratio of damping factor for the finite 
disc to that for the infinite disc, equation JlOt , is roughly given by 
a phenomenological factor. 



(17) 



I^M('^)linflmtc 

This also provides a good fit to the numerical results [though worse 
than that using the Fourier transform of equation J16H, as shown in 
Fig. [3] 

In order to explain the rather significant difference between the 
form of the Green's function between the finite and infinite disc, 
power law vs. exponential, we have solved numerically the case 
of mass supply at an intermediate radius, ro <c ro„t, where 

''out tS 

the outer radius of a finite disc. We have assumed rom/rin = lO"" 
and n = 3/4. We have found that the exponential behaviour of 
the response to a pulse is due to the finite size of the disc. The 
Green's function first follows the form of that for the infinite disc, 
with a power law decay, but it starts to depart from the infinite- 
disc solution at ?/fvisc(''o) ~ I'out/^o and then it decays exponentially. 
This is illustrated in Fig.|4] for two choices of rom/fo- The-power 
law asymptotic form is due to the assumption of the infinite r^^t, 
allowing the mass flow to diffuse back to the inner radius back from 
arbitrarily large radii. 

We also note that the overall form of the finite-disc Green's 
function, not only its exponential decay, appears common for dif- 
fusion problems. For example, the Green's function for Thomson- 
scattering diffusion of photons out of a spherical electron cloud has 
a very similar shape to that found for the present problem, see, e.g., 
IZdziarski^isra & Gierlinskil j2009l) . 



3 FURTHER NUMERICAL MODELS 

We present most of numerical results in this section for n = 3/4, 
appropriate for outer disc regions, where most of the viscous damp- 
ing takes place. We assume Mout = Mm at f = 0, and we start to 
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Figure 5. Time evolution of tlie mass accretion rate for a step-function in- 
crease of tlie mass supply rate at the disc outer edge from Mom = up to 1 
at ; = 0, for ^oui/ri,, = 10^-10^. These results are compared to those for the 
infinite disc, equation )13t . shown by the thick solid curve. 
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Figure 6. Time evolution of the mass accretion rate at fi = 3/4 and r„ut I I'm = 
K)-* for the case of the sinusoidal modulation with A(rout) = 1-8, shown by 
the dotted curve, turned on at ? = 0. The solid and dashed curves show 
M\n(t) for ?visc('"oui) = P and 0.4P, respectively. 



change Mout at / > 0. We consider first the mass supply rate in- 
creasing as a step function from M^ut = to 1. Fig. [5] shows the 
resulting evolution of the mass accretion rate. The inner mass ac- 
cretion rate increases on the time scale of ~ ?visc(''out)- The steady 
state is reached after a few /visc(''out)- The accretion rate reaches an 
asymptotic temporal profile independent of rout/Hn at its very large 
values, which is due to most of the viscous damping taking place 
in outer parts of the disc. We find the steady state is achieved much 
faster than in the case of an infinite disc, which is due to its Green's 
function being much broader than the one for finite disc, see Figs. 

Then, we consider periodic modulation of the mass supply 
rate. After turning on the modulation, the disc settles down into 
a steady-state periodic behaviour after some number of fvisc(''out)- 
We define the steady-state fractional variability amplitude. 



A(r): 



(M) 



(18) 



where Mmax('') and Mn,in(r) are the maximum and minimum flow 
rates, respectively, and (M) is the average rate. 
We first consider a sinusoidal modulation, 



(M) 



A(ro„t) . 2nt 

1 + sin — 

2 P 



(19) 



where P is the modulation period (usually, but not necessarily, 
equal to the binary period, Poit)- Fig. [6] shows an example of the 
disc settling into a periodic behaviour after starting the modulation 
at f = 0. The inner mass accretion rate shows regular sinusoidal 
behaviour after several ?visc(''out)- Fig. |7] shows the fractional vari- 
ability as a function of the radius for this case for several values 
of t^hc(roui)/P- We see that the final value of the modulation am- 
plitude, Ain, is reached already at r » r,„. Obviously, the viscous 
damping becomes stronger with the increasing /visc(''out)/^'. 

In the case of a large eccentricity, accretion can take place 
only around periastron, as it is the case, e.g., in Be X-ray binaries. 
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Figure 7. The amplitude of the steady-state periodic variability of the 
mass flow rate as a function of the radius, A(r), for a number of values of 
'visc(''oui)/P- The parameters are n = 3/4, rout/''iii = 10'* and the amplitude 
of the sinusoidal modulation is A(roui) = 1.8. 



We thus approximate the disc supply rate as a periodic rectangular 
function. 



M, 



:<M) 



P/AP, k<t/P <k + AP/P; 



(20) 



0, otherwise, 
where AP is the duration of the accretion, and k is an integer. We 
also consider Mam corresponding to a Roche-lobe overflow, 



27rf , 

M„„, = Ml exp ( C sin — | , 



(21) 



see, e.g., IZdziarskietd] ( 12007 ah . who obtained C - 1.5 in the case 
of the ultracompact low-mass X-ray binary 4U 1820-303 at the 
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Figure 8. The amplitude of the periodic variability of the mass accretion 
rate at the inner edge of the disc for AP/P = 0.01 (thick solid curve) and 
0.1(thin solid curve), for C = 10 (thick dashes) and 1 (thin dashes), and 
for sinusoidal case with A(rom) = 1.8 (dots). In all cases, n = 3/4 and 

'-out/nn = 10"*. 



maximum of its modelled eccentricity. Here Mj is a normalization 
factor. Fig. [8] shows the relative variability amplitude at the inner 
edge as a function of ?vi.sc(''ou()/^' for sinusoidal, periodic rectangu- 
lar and Roche-lobe overflow cases. When ?visc(''oiit) ^ 10 P, the am- 
plitude is very strongly damped and A(r,„) < 0.01 even if P„yt,/AP 
or C are very large. The form of the damping is qualitatively similar 
in all cases. Thus, our result that the amplitude is strongly damped 
for a large /visc(''out)//' does not depend on the functional form of 
the mass supply rate. 



4 APPLICATIONS TO X-RAY BINARIES 

The viscous time is, in the a prescription, given by 

fvi»c = -^-f— 1 , (22) 
aLlK \ vkI 

where vk is the Keplerian velocity (and Cs/vk = H/ r, where H is the 
disc scale height). The period of the variability of the mass supply 
rate, P, usually equals the orbital period, which is given by the Ke- 
pler law, Porb = 2n ^la^ IGMi{\ + q), where q = MijMu and Mi, 
M2 and a is the mass of the compact star, the donor mass, and the 
semi- major axis, respectively. Then (see also iGilfanov & ArefievI 
l2005h . 



('"out) (1+?) 



1/2 



Vk/ \ a I 



(23) 



In the standard thin accretion disc, the ratio Cs/vk is very small, 
~ 0.02. In low-mass X-ray binaries, the size of the disc is compa- 
rable to the size of the Roche lobe around the compac t object, Ri 
jPapaloizou & Pringlelll977l : lGilfanov & Arefievll 20051) , which, for 
e <K \ and g' ^ 1 characteristic to that class of system, is in turn 
comparable to a. Thus, rom/a ~ 0.5, and then the viscous time is 
much longer than the orbital period, fvisc(''ou()/^'orb ^ 10^/0-. Then, 
our results imply that any periodic variability of the mass accre- 
tion rate at the inner disc is completely negligible in the standard 



model even if the mass supply to the outer disc is periodic. This 
implies, in particular, that the ~ 1 per cent orbital modulation seen 
in X-rays in 4U 1820-303 is not due to periodic ity of the mass 
supply predicted by the triple model of that system l lZdziarski et al.l 
l2007al) , provi ded the accretion disc is standard. This confirms the 
supposition of Zdziarski et al ]| |2007bl) . On the other hand, the vis- 
cous damping of the superorbital variability in 4U 1820-303, with 
P 170 d, is negligible unless a ^0.01. 

We also note that our results rule out the eccentric low-mass 
X-ray binary model of Maccarone (2 00j) f or the X-ray flares 
in NGC 7469 found by Sivako lf et all ( 120051) . iMaccaronel ( |2005|) 
claims that since the flares are super-Eddington, the accretion discs 
are radiation-pressure dominated and then ~ vk and t^\^^ < P„^f,. 
However, super-Eddington discs are radiation-pressure dominated 
only in inner regions whereas the outer parts are still gas-pressure 
dominated and geometrically thin (unless the accretio n rate is ex- 
tremel y high). Specifically, for parameters suggested bv lMaccarond 
l l2005h . Ml = 1.4Mo, Ml = IMq, P„b = 15 hr, a = 0.3, the disc 
size of 0.9 of the Roche lobe size, Mc-/Le = 100 (where is 
the Eddington luminosity), we have obtained t^isc ~ 10" Porb using 
standard disc equations. Obviously, there will be then no observ- 
able enhancement of th e X-ray flux ar o und th e periastron. 

We also note that ISivakoff et all ( l2009h have found no peri- 
odicities in X-ray flares in NGC 7469, which provides independent 
evidence against the eccentric binary model. The flares may instead 
have similar nature to the aperiodic flares discover ed from the black 
hole b inary Cyg X- 1 with a range of time scales jGies et al ] |2003l; 
iGierli iiski & Zdziarski 2003). 

"Gilf anov & ArefievI ( 120051) found breaks in the power spec- 
tra in a number of low-mass X-ray binaries at frequencies propor- 
tional to the orbital frequencies. They interpreted those frequencies 
as equal to Fj'^, which then implied that t^\^^ were much shorter than 
that predicted by the standard model, which effect they explained 
by the presence of hot coronal flows covering outer parts of the 
disc. In particular, they claimed t^^^^ ~ 300 s ~ Porb/2 in 4U 1820- 
303. After such a short time interval, most of the large, > 10, mass- 
supply modulation at the outer disc due to eccentricity predicted by 
the triple model (Zdziarski et al. 2007a) would still be present close 
to the neutron star surface. However , such modulation is clearly not 
seen. Thus, t^\^c ~ 300 s claimed bv IGilfanov & ArefievI ( l2005i) for 
4U 1820-303 is ruled out in the framework of this model. On the 
other hand, it may be that the cold disc is covered by a fast coro- 
nal flow having a short viscous time and carrying ~1 per cent of 
the mass flow, which could then account for the observed X-ray 
periodicity. 

Then, high-mass X-ray binaries have usually ^ » 1 , for which 
the size of the Roche lobe of the compact object in the circular case 
is (Paczyiiski 1971.) . 



(24) 



1/3 • 



If 

''out were ~ R\ (as expected for Roche lobe overflow), the time 
scale ratio of equation \22>\ would become oc (f, and the conclu- 
sion of fvisc( ''out) /Porb » 1 would Still hold. Taking into account 
eccentricity can be done in a simplified way by a substitution of 
a — > a(l - e) (i.e., using the periastron separation as constraining 
the disc size), which reduces the time scale ratio by (1-e)''-, which 
is significant only for e ~ I. 

However, discs in wind-accreting systems may be truncated at 
radii much smaller than the Roche lobe radii of the compact ob- 
ject, e.g., Shapiro & Lightman (1976). The outer edge of the disc 
is then determined by the specific angular momentum in a Keple- 
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rian orbit around the compact object, (GMifou,)''^, being equal to 
the angular momentum carried by the accreting gas, (l/2)Voibrl/a, 
where r.^ ^ 2GMi/vJ^j is the accretion radius, Vo,.(, is the relative 
velocity of the stars, Vrei is the relative velocity of the compact ob- 
ject and the vyind, vL = vl ^, + v^jj,^, and Vwmd is the wind velocity 
dShapiro & Lightmanll 19761) . This yields, 

= ^Jl (25) 

a b\\+qf{2a-df' 

where d is the orbital-phase dependent distance between the stars 
and b = Vrei/Vo rb- Notc that equa tion V.5\ can be used only if it 
yields < /?, . IPravdo & GhoshI JioOl) assumed that the neutron 
star in GX 301-2 accretes from a slow, dense, circumstellar wind, 
in which case b ^ \. Then at periastron, ronilci becomes 4[(1 - 
e)/(l + q){\ + e)]\ which for ^ ^ 20, e ^ 0.5 (Kaper at al. 2006) is 
2 X 10"^ Equation then yields /visc(''out)/^'mb ~ 2 x 10"'*/q'. 
Thus, the viscosity damping of the orbital periodicity is negligible 
in GX 301-2. Similar considerations apply to eccentric Be/X-ray 
binaries. 

The case of Cir X-1 is less clea r given the uncertainty about 
its system parameters. IClarkson et alj j2004h have modelled the or- 
bital phase d ependence of the mass transfer rate in Cir X-1 using 
the model of iBrown & Bovlj ( Il984h . and they estimated /vise ~ 10 
d, i.e., also less than its Po,b- Thus, viscosity has probably a minor 
effect on the observed orbital periodicity in that object. Observa- 
tionally, we do see strong periodicity at Po,!,, which also implies 

We also note that the time scale of state transitions in black- 
hole or neutron-star binaries is likely to be the viscous time scale at 
the radius of the mass supply (different from the case of outbursts 
of either X-ray o r dwarf novae, where the time scale is by a factor 
of ~ HI r shorter, lMeveilll984 iKato et al.|[T993) . Then, our results 
appear to rule out the mass supply in Cyg X-1 (with Porb = 5.6 
d) being at the outer edge of a standard disc with the size compa- 
rable to its Roche lobe radius. If this were the case, fyis would be 
~ lO'/'orb/o'. Then, any change of the mass supply rate to the disc 
would result in a change of the mass accretion rate in the inner disc 
being on a time scale of » Paib- This is much longer than the ob- 
served time s cale of state transitions in th is system of a few PoA 
(e.g., fig. 5 in lZdziarski & Gierlinskill2004 'l Thus, the state transi- 
tions in Cyg X-1 cannot be caused by a change of M at the outer 
edge of a standard accretion disc. Indeed, though the donor almost 
fills i ts Roche lobe, the wind accretion is still dominant ( Gies et al. 
l2003h . and the size of the disc is then likely to be much smaller 
than that of the Roc he lobe of the b lack hole. Equation l l25b with 
q^2.?,± 0.4, b~2 dCies et alj2003l) and e = yields then a rather 
small disc, with r„^ila ^ 3 x 10""*, at which equation ([23} yields 
fvisc(''out)/^oib ~ 0.004/0-. (Note that the wind velocity around the 
compact object in Cyg X-1 is not well determined, which makes the 
determination of ro„t/a rather uncertain.) The combination of the 
actual disc size in Cyg X-1 (which may be intermediate between 
those of the pure Roche-lobe overflow and wind accretion cases) 
and the viscosity parameter may correspond to fvisc(''out) ~ ^orb> as 
it appears to be implied by the observations. Still, a large accretion 
disc can be present if the actual viscous t ime is much shorter than 
the standard one, as in the model of Gilfa nov & ArefievI j2005l) . 

We can also invert equation ( 123b in order to obtain the disc 
size corresponding to a given /visc(''out)/^'orb. This yields, 

(26) 



{2naf^^ 



(1+9) 



1/3 



(If 


^visc(''out) 


2/3 




P orb 





We see that fvisc(''oiit) ~ ^"0* implies rom/a 1 in general. 



Throughout the present study, we have assumed that a disc 
has an axisymmetric structure and ha ve examined only the v iscous 
damping effect. On the other hand, iHavasaki & Okazakil j2005l) 
performed three-dimensional smoothed particle hydrodynamics 
simulations of Be/X-ray systems and claimed that time-dependent 
mass transfer to the disc will produce a non-axisymmetric struc- 
ture because of the ram pressure by the incoming stream. Then, 
a spiral wave is created and its inward propagation significantly 
enh ances the mass-accretion rate onto the neutron star (see also 
Hava saki. Mineshige & Hoi |2008| in the context of binary black 
holes). If this is the case, it is possible to produce appreciable light 
variations even when the viscous time scale at the outer edge of 
the disc is much longer than the time scale of the mass input. Our 
simple analyses presented in this paper does not allow to exam- 
ine this effect and multi-dimensional simulations are necessary in 
future work. 

In our applications, we have considered only binaries. How- 
ever, our results on the viscous damping of variability hold for any 
accretion disc, also those around supermassive black holes in active 
galactic nuclei. 



5 CONCLUSIONS 

We have calculated in detail the effect of viscous damping of vari- 
ability in accretion discs. We have first considered the infinite-disc 
model of Lynden-Bell & Pringle ( 1974). For it, we have obtained a 
simple form of the Green's function for M, equation Q, as well as 
its Fourier transform, equation l llOt . The latter gives the fractional 
viscous damping at a given variability frequency, approximately 
given by ~exp(- ^fo)). As a simple application of our Green's func- 
tion, the case of a step function increase of M yields an approxi- 
mately exponential approach to the steady state, ~exp(-l/T). 

We have then calculated numerically a realistic finite-disc 
model with the mass supply at its outer boundary. We have found 
significant differences with respect to the infinite-disc model. In 
particular, a single pulse at the disc outer edge results in an expo- 
nential decay after a fraction of the viscous time at the disc outer 
edge, whereas the infinite-disc model yields a much broader power- 
law decay. The power-law behaviour is due to the infinite size, al- 
lowing the flow to diffuse back from arbitrarily large radii. Given 
the differences, we have fitted the finite-disc Green's functions by 
a simple form, and calculated its Fourier transform. We have also 
calculated numerically time-dependent disc behaviour for a num- 
ber of more complex cases, in particular for periodic modulations 
of the mass supply rate. 

We have applied our results to a number of X-ray binaries. Pe- 
riodic variability resulting from eccentricity is usually very strongly 
damped in low-mass X-ray binaries, which accretion discs have the 
size comparable to that of the Roche lobe of the compact object. 
This is the case, in particular, in 4U 1 820-303 or in flaring binaries 
in NGC 7469. 

The periodic mass transfer in eccentric binaries is much easier 
to observe in high-mass X-ray binaries, where accretion is via a 
wind and the accretion discs can be truncated at radii much smaller 
than the size of the Roche lobe, resulting in the viscous time being 
relatively short. In particular, the viscous damping is negligible in 
GX 301-2, Be/X-ray binaries and Cir X-1. In the case of Cyg X- 
I, the observed time scale at which state transitions take place is 
comparable to its orbital period, which implies the accretion disc 
being also truncated at a radius much lower than the Roche-lobe 
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size. This may be due to the presence of wind accretion in this 
system. 
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APPENDIX A: VISCOUS PROPAGATION OF M 
BETWEEN TWO DISC RADII 

The explicit form of G^Cr, t) of equation (O for yu > is, 

2yu2(^"'' + 1) 



Gm(^,t) = exp 



X I, 



.2 £ 1/(2^) 



r 

2tl/(2,i) 



At n = \ ill = 1 12), 



whereas its Fourier transform is. 



|cosh - - ^sinh -j , 



(Al) 



(A2) 



(A3) 
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